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Abstract 

We study a catalytic branching process (CBP) with any finite set of catalysts. This 
model describes a system of particles where the movement is governed by a Markov chain 
with arbitrary finite or countable state space and the branching may only occur at the 
points of catalysis. The results obtained generalize and strengthen those known in cases 
of CBP with a single catalyst and of branching random walk on Z'^, d G N, with a finite 
number of sources of particles generation. We propose to classify CBP with N catalysts as 
supercritical, critical or subcritical according to the value of the Perron root of a specified 
N X N matrix. Such classification agrees with the moment analysis performed here for 
local and total particles numbers. By introducing the criticality set C we also consider 
the influence of catalysts parameters on the process behavior. The proof is based on 
construction of auxiliary multi-type Bellman-Harris processes with the help of hitting 
times under taboo and on application of multidimensional renewal theorems. 

Keywords and phrases: catalytic branching process, classification, hitting times under 
taboo, moment analysis, multi-type Bellman-Harris process. 
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1 Introduction 

We consider the model of catalytic branching process (CBP) with a finite number of catalysts. It 
describes a system of particles moving in space and branching only in the presence of catalysts. 
More exactly, let at the initial time t = 0 there be a single particle that moves on some finite 
or countable set S according to a continuous-time Markov chain r] = {r]{t),t > 0} generated 
by infinitesimal matrix Q = {q{x,y))x,y£S- When this particle hits a finite set of catalysts 
W = {tci,..., wn} C S, say at the site Wk, it spends there random time having the exponential 
distribution with parameter fik > 0. Afterwards the particle either branches or leaves the site 
Wk with probabilities and 1 — ak (0 < < 1), respectively. If the particle branches (at the 

site Wk), it dies and just before the death produces a random non-negative integer number ^k of 
offsprings located at the same site Wk- If the particle leaves Wk, it jumps to the site y Wk with 
probability —(1 — ak)q{wk,y)q{wk,Wk)~^ and continues its movement governed by the Markov 
chain rj. All newly born particles are supposed to behave as independent copies of their parent. 

We assume that the Markov chain rj is irreducible and the matrix Q is conservative (i.e. 

~ 0 where q{x,y) > Q ioi x y and q{x,x) G (— cxd,0) for any x G S). Denote 
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by fk{s) := Es^*, s G [0,1], the probability generating function of k = 1,... ,N. We will 
employ the standard assumption of existence of a hnite derivative /((I), that is the hniteness 
of for any k = 1,..., N. 

Study of CBP with a single catalyst was initiated even in the XX century (see, e.g., 0). In 
this regard we also mention a recent paper na where the main tool for the moment analysis 
of the process was the spine technique, i.e. “many-to-few lemma”, and renewal theory. Note 
that the generalization to an arbitrary finite set of catalysts is not straightforward, since there 
exists a ’’competition” between catalysts. An important special case of several catalysts (where 
S = d G N, and the Markov chain r] is symmetric space-homogeneous random walk 
with hnite variance of jump sizes) was examined in [33]. There the sufficient conditions for 
exponential growth of the particles numbers were obtained by analyzing spectral properties of 
evolution operators. We propose another way permitting to establish more precise results. 

For branching processes the natural interesting problem is the analysis of asymptotic be¬ 
havior (as t —)• oo) of the local and total size of population at time t (for various branching 
processes without catalysts see, e.g., 1271 ). Let /i(f) stand for the total number of particles 
existing in CBP at time t > 0. In a similar way we dehne local numbers /x(f; y) as quantities of 
particles located at separate points y E S at time t. In this paper our aim is three-fold. Firstly, 
we introduce a classihcation of CBP (with N catalysts) treating it as supercritical, critical or 
subcritical whenever one has, for the Perron root p of a certain N x N matrix, p > 1, p = 1 or 
p < 1, respectively. Moreover, the criticality set C revealing the inhuence of catalysts strength 
is introduced and characterized. Secondly, we implement the moment analysis of the local and 
total particles numbers to justify the naturalness of the proposed classihcation (indeed, one 
will see that the asymptotic behavior of the moments of any order is determined essentially by 
the introduced class of CBP). Thirdly, we consider some particular cases of CBP and thereby 
show that our study not only generalizes the results in previous works but even rehnes them. 

Our approach consists in involving hitting times under taboo (see, e.g., 0. 0 and HU, 
Ch.2, Sec.11) and construction of auxiliary Bellman-Harris branching processes with at most 
A^(A^-|- 1) -l- 1 types of particles. This approach is inspired by [HI| where the branching random 
walk on Z with a single catalyst was investigated by means of introducing hitting times (without 
taboo) and a due two-type Bellman-Harris process. Note also that in [101, for the study of a 
discrete-time branching random walk on Z with multiple catalysts, the authors considered an 
embedded multi-type Galton-Watson branching process, resulting in ’’forgetting/erasing the 
time spent between catalysts”. The latter approach is fruitful for classihcation of branching 
random walk since multi-type Galton-Watson and multi-type Bellman-Harris processes have the 
same supercritical, critical or subcritical regimes. However, for subsequent study of branching 
random walks (moment analysis, limit theorems etc.) multi-type Galton-Watson process is not 
sufficient and thus Ph.Carmona and Y.Hu used in [10] another technique such as ’’many-to- 
few lemmas”. Our hrst auxiliary process is constructed in such a way that the study of CBP 
can be mainly reduced to analysis of the Bellman-Harris process. The number of particles in 
this process is chosen to guarantee its indecomposability and cannot be less than N and need 
not be greater than N{N -|- 1). It is well-known (see, e.g., [27], Ch.4, Sec.5, 6 and 7) that 
an indecomposable multi-type Bellman-Harris process is classihed as supercritical, critical or 
subcritical according to the value of the Perron root of the mean matrix. This is the foundation 
for classihcation of CBP. Furthermore, LemmajUbelow gives even more convenient classihcation 
with the help of the Perron root of a specihed irreducible N x N matrix. One more auxiliary 
Bellman-Harris process is indispensable for study of total particles numbers in CBP in case 
of transient Markov chain r]. This process is taken to be decomposable with a hnal type of 


2 


particles. The treatment of the auxiliary Bellman-Harris processes allows us not only to classify 
CBP but also to derive a system of renewal equations for the mean local and total particles 
numbers in CBP. Afterwards, to implement the moment analysis of the local and total particles 
numbers we use multidimensional renewal theorems established in na and [23] . 

This approach has advantages since there is elaborated theory of multi-type Bellman-Harris 
processes and a vast majority of its results can be applied to the auxiliary Bellman-Harris 
processes leading to the new results for CBP. Mention in passing some recent works on multi¬ 
type Bellman-Harris processes, see, e.g., [IB], (30] and [32] . 

Observe also that CBP can be considered as a Markov branching process with at most 
countably many types of particles, since the location of a particle can be associated with its 
type. Theory of branching processes with countably many types of particles, despite of its 
long history (see, e.g., [23|), has not been systematized until now in view of its complexity. 
Some new papers such as 0. @1. nn] and [22] make an important contribution to this research 
direction. However, the study there does not cover the results presented in our paper. Among 
other investigations of models describing particles movement and breeding we refer to recent 
papers 0. mi and [21]. Branching random walks analyzed there are homogeneous in space 
whereas the main feature of CBP is spatial non-homogeneity. Some differences in behavior 
of a homogeneous branching random walk on Z and its catalytic counterpart are discussed in 
[In] . Concluding the introduction we mention the close relation between a catalytic branching 
random walk on Z and a super-Brownian motion with a single point catalyst (see, e.g., na 
and [T5]h 

Now we describe the structure of the paper. In section we specify the hrst auxiliary 
Bellman-Harris process and propose a classihcation for CBP. Section is devoted to solution 
of the following problem. Assume that we £x the Markov chain generator Q in CBP and vary 
“intensities” of catalysts, e.g., nik ■= k = 1,. .., A^. What is the set C C R+ such that 
CBP is critical iff m = (mi,... ,mjv) € Cl In particular, what is the proportion of “weak” and 
“powerful” catalysts (possessing small or large values of m*, respectively) in critical regime? In 
section we obtain a complete description of C by means of equation involving determinants 
of some matrices, indicate the smallest parallelepiped [0, Mi] x ... x [0, M^r] containing C 
and illustrate our approach by two pictures of the set C when N = 2 and N = 3. Observe 
that these plots are not a product of simulation and display the result of direct computation. 
Section contains Theorem and its proof representing the moment analysis of the local 
and total particles numbers in CBP. There we construct the second auxiliary Bellman-Harris 
process required for study of the total size of population whenever Markov chain rj is transient. 
Section demonstrates applications of our results to the catalytic branching random walk on 
Z*^ and to the branching process with a single catalyst. A detailed comparison of our results 
with those known before is given as well. 

2 Auxiliary Bellman-Harris process 

Let us briefly describe a Bellman-Harris branching process with particles of n types, n G N. It 
is initiated by a single particle of type i = 1,..., n. The parent particle has a random life-length 
with a cumulative distribution function (c.d.f.) Gi(t), t > 0. When dying the particle produces 
offsprings according to a probability generating function gi{s), s = (si,..., Sn) G [0,1]"'. The 
new particles of type j = 1,... ,n evolve independently with the life-length distribution Gj{-) 
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and an offspring generating fnnction gj{-). Let 

M ■= with rriij = 5., c/4=(i,...,i) 

be the mean matrix of the process. The Bellman-Harris branching process is called indecom¬ 
posable if the non-negative matrix M is irreducible (for the latter notion, see, e.g., [26], Ch.l, 
Sec.3). Moreover, if its Perron root p (i.e. eigenvalne having the maximal modulus) is such 
that /9>l,p=lorp<l, then the Bellman-Harris process is called supercritical, critical or 
subcritical, respectively (see, e.g., [27|, Ch.4, Sec.5, 6 and 7). Denote the number of particles 
of type j existing at time t by t > 0, j = 1 ,..., n. 

Before demonstrating how an auxiliary Bellman-Harris process can be constructed in the 
framework of CBP we have to introduce some notation. Consider a particle moving on the set 
S in accordance with the Markov chain generated by Q and starting at state x. Let H^x,y, 
x,y E S, H (Z S, he the time spent by the particle after leaving the starting point x until the 
hrst hitting y if the particle’s trajectory does not pass H. Otherwise (if the particle’s trajectory 
passes H before the hrst hitting ?/), HTx,y = oo. The (extended) random variable HTx,y is called 
a hitting time of state y under taboo H after the hrst exit out of the starting state x. Denote by 
H^^^yit), t > 0, the improper c.d.f. of HTx,y Evidently, HFx,y{0) = (4,y - ^)q{x,y)q{x,x)~^ 
where is the Kronecker delta. Explicit formulae for the probability of hniteness of H^x,y, i-e. 
for HFx,y{oo) = limt_^oo HFx,yit), via taboo probabilities and Green’s function were derived in 
[9]. Whenever the taboo set H is empty we write Tx,y and Fx,y{-) instead of (nfx^y and (/iFx,y{-), 
respectively. 

Return to CBP. We tentatively assume that CBP starts at Wp for some p = 1,..., N. Set 
Kj := {k = 1,..., N : w^F^.^^^ioo) - > 0} for Wk:=W\ {wk} and j = 1,..., N. 

Let Kj = {k{l,j),...,k{\Kj\,j)} where 1 < k{l,j) < k{2,j) < ... < k{\Kj\,j) < iV and | ■ | 
is the cardinality of a hnite set. We divide the particles population existing at time t > 0 into 
L-\-l groups with L := + \Fj\- The particles located at time t at Wj form the j-th group 

having cardinality p{t]Wj), j = 1,N. Consider a family consisting of particles which have 
left Wj at least once within time interval [0,t], upon the last leaving Wj have not yet reached 
W by time t but eventually will hit Wk(ij) before possible hitting H4(ij), i = ■ ,\Fj\, 

j = 1,...,A^. This family has cardinality denoted by Pjj{t) and corresponds to the group 
number (L(j) -|- i) where L[j) ■= N -\- The group number L -\- 1 comprises the rest 

of particles not included into the above L groups. Note that the last group consists of the 
particles having inhnite life-lengths since after time t they will not hit the set of catalysts any 
more. So, after time t these particles will not produce any offsprings and have no influence on 
the numbers of particles in other L groups. 

Now we can introduce an auxiliary Bellman-Harris process to employ it for the study of 
CBP. Consider an L-dimensional Bellman-Harris process starting with a single particle of type 
p and having the following c.d.f. and offspring generating functions 


N 


Fljit') 1 e , gj{^') Q^j) ^ w^FWj 


k=l 

\Kj\ / N 

-|-(1 — Tjj(cx))s2,(j)+i -|- (1 — ttj) ( 1 ■“ ^^WkFwj,w^.{oo) j , 

i=i \ k=i 

Tij(t) 


( 1 ) 


( 2 ) 
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where Tij{t) := 0 < t < oo and z = l,...,\Kjl 

j = l,...,iV. It is easy to see that for the process constructed in this way the vector 
{Ziit),..., Ziit)), at each t > 0, has the same distribution as the vector whose j-th com¬ 
ponent is /i(t; Wj) and (-h(j) + z)-th component is i = 1, ..., \Kj\, j = 1, ..., N. 

Mean matrix M = of the introduced Bellman-Harris process is of block form, i.e. 


M = 


( Mi^i Mi^2 

^ 2,1 ^ 2,2 

V Mtv+I,! Mn+1,2 


Ml, 


N+1 


\ 


M 


2,7V+1 


(3) 


d^V+l,Af-|-l / 


Here the A^xA^ matrix Mi,i = {jnk^i)ki=i possesses the following entries rrik^i = + 

otk)wiFWk^wiiS^)■ For i = 1,..., A^, the elements of A^ x \Kj\ matrix Mij+i vanish everywhere 
except for the j-th row and for z = 1,..., \Kj\ one has = (1 — aj)Tij{oo). When 

z = 1,..., A^ the lATjl x N matrix Mj+i,i is such that omitting its j-th null columns for all 
j ^ Ki one gets the \Ki\ x lATjl identity matrix. For i,j = 1,... ,N, the \Ki\ x \Kj\ matrices 
Mj+ij+i have zero entries. 

Let us verify that the specihed Bellman-Harris process is indecomposable. In view of ir- 
reducibility of the Markov chain rj there exists a finite path from Wk to wi with a positive 
probability for each k,l = 1,..., N. In the framework of CBP such path has also positive prob¬ 
ability, since being at catalyst site Wi a particle can leave it without branching with positive 
probability I — ai for each z = 1,..., A^. Among the sites visited successively by that path 
of T] choose those from W, say, znp(o), znp(i),..., znp(j), znp(j+i) with p(0) := k and p( J -|- 1) := 1. 
By the construction of auxiliary Bellman-Harris process this path corresponds to transforma¬ 
tions of a particle type from k to 1. Namely, if the path of r] hits znp(j+i) immediately after 
leaving znp(j) then the particle type in Bellman-Harris process changes from p{j) to p{j -|- 1). 
Otherwise, p{j -|- 1) G Kp^p i.e. p{j -|- 1) = k{i,p{j)) for some z = 1,..., |A'p(j)|, and the 
particle type change from p{j) to p{j + 1) involves the intermediate type L{p{j)) + i. Hence, 
we conclude that, for each k and I from {!,..., A^}, there exists n = n{k,l) G N and a col¬ 
lection {r(l),..., r(rz)} C {1,..., L} such that mfc,r(i)^r(i),r( 2 ) • • • > 0. Since by matrix 

M dehnition in (|^ one has rnL(^j)^i^k{i,j) = 1 > 0 and mj^L{j)+i > 0 for each z = 1,..., \Kj\ 
and j = 1,..., A^, the previous statement holds true for /c, Z = 1,..., L as well. Consequently, 
we have checked that for each k and I from {!,..., L}, there exists n = n{k, /) G N such that 
m^i^i > 0 where is the {k, /)-th element of M'^. So, we have verihed that the matrix M is 
irreducible. Moreover, if we consider the nontrivial case (1) > 0 for some z = 1,..., A^, then 
rrii^i > 0 and the matrix M is acyclic (or aperiodic) (see, e.g., |26], Ch.l, Sec.2). Thereby, in 
this case we have shown that M is a primitive matrix (see, e.g., [2E], Ch.l, Sec.3). 

On account of the Perron-Frobenius theorem for irreducible matrices (see, e.g., [2S], Ch.l, 
Sec.4) the matrix M has a positive real eigenvalue p{M) of maximal modulus which is called 
the Perron root of M. Using the classification of the constructed Bellman-Harris process we 
call CBP supercritical, critical or suhcritical if p(M) > 1, p(M) = 1 or p(M) < 1, respectively, 
where M is specified by (|^. Since M is independent of the CBP starting point, the same is 
true for the classification of CBP. So, from here on we omit our assumption that CBP starts 
at some zUp G lU. 

Let G*(A) := e~^^ dGi{t), A > 0, be the Laplace-Stieltjes transform of c.d.f. Gi and 

H{X) := {G*{\)mij)fj^^. Note that H{0) = M. Matrix H{X) is irreducible in view of irre- 
ducibility of M. Put also D{X) = {dij{X))F^i, A > 0, with dij{X) := aif-{l)G*{X) -|- (1 — 
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ai)G*{X)wjF*^.^wj{^)- It is easily seen that the matrix D = {dij)F^i dehned as D := -D(O) has 
entries 

(I Oli)WjF 11^,111 jis^') ■ (4) 

For an irreducible matrix A, let p{A) stand for the Perron root of A. The following statement 
gives a convenient criticality condition for CBP. 

Lemma 1 For each A > 0, the matrix -D(A) is irreducible and the values p{D{X)) and p{H{X)) 
are simultaneously greater than 1, equal to 1 or less than 1. In particular, the Perron roots 
p{D) and p{M) are either both greater than 1, or equal to 1, or are less than 1. 

Proof. The irreducibility of D (and, therefore, of D{X)) is established in the same manner 
as for matrix M. In accordance with the Perron-Frobenius theorem matrix H{X) has a strictly 
positive left eigenvector v corresponding to the eigenvalue p{H{X)). Consequently, H{X)~^\ = 
p{H{X))\ where T means transposition and v is considered as a column vector. In other words, 
(/f(A)^ — p{H{X))l) V = 0 for / being the identity matrix and 0 standing for the zero of R^. 
Apply some equivalent transformations to the obtained system of equations, namely, for each 
j, k = 1,..., N add the (T(j) + i)-th row multiplied by p{H{X))~^ to the k-th row if 

k G Kj, i.e. k = k{i,j) for some i = 1,..., \Kj\- Focusing on the first N equations only, we 
deduce that 

vj {Hu(\) + p(H(X))-'R(X) - p(H(\))I) = or 

Here we denote Hi^i{X) := {G*{X)mij)F^^ and matrix R{X) := {rij{X))F^^ has the following 

entries rjj(A) := (1 - ai)G*{X) (^WjF*wi,wjW “ whereas vq is a strictly positive 

vector consisting of the first N coordinates of v and 0 means the zero of R'^. Note that 
D{X) = Hi^i{X) + R{X). Evidently, if p = p{H{X)) > 1 then 

(D(A) - pU) = vj + p-'B(A) - pi) 

< vj (i/,,,(A)+p-‘fl(A)-p/) < vj (Hu(\)+R(\)-pI) = vj (D(A)-p/). 

In a similar way, if p{H{X)) < 1 then the previous chain of inequalities holds true when each 
symbol < is replaced by >. Finally, we get 

^(A)"^vo > p(iP(A))vo and T>(A)’^vo < p{H{X))\q for p{H{X)) > 1, (5) 

D{X)~^'Vf) < p(iP(A))vo and iA(A)'''vo > p(i/(A))^vo for p{H{X)) < 1. (6) 

On account of Theorem 1.6 in [2S], Ch.l, Sec.4, the second relation in ([^ and the first one in 0 

imply that p{D{X)) < p{H{X))'^ for p{H{X)) > 1 and p{D{X)) < p{H{X)) for p{H{X)) < 1. Ex¬ 
amining the proof of the Perron-Frobenius theorem we conclude that due to the first inequality 
in ([^ and the second in (|^ one has p{D(X)) > p{H(X)) for p(H{X)) > 1 and correspondingly 
p{D{X)) > p{H(X))^ for p(H{X)) < 1. Consequently, if p{H(X)) > 1 then p(H(X)) < p{D{X)) < 
p(i/(A))^ whereas if p{H{X)) < 1 we obtain p{H{X)Y < p{D{X)) < p{H{X)). These estimates 
entail the assertion of Lemma [U □ 

Definition 1 CBP is called supercritical, critical or subcritical whenever p{D) > 1, p{D) = 1 
or p{D) < 1, respectively, where D is the matrix specified by way of (|^. 
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3 Structure of the criticality set 

In this section we focus on the nonsingular case G (0,1) for each i = One 

can see that, given a Markov chain generator Q and a collection ai,... ,aN, the matrix D 
depends on = //(I) only and not on the explicit form of i = 1 ,..., At, that is 

D = D{mi,... ,mjv). Dehne the criticality set as 

O := {(mi,..., miv) G R+ : p{D) = 1}. 

Study of the criticality set structure allows us to understand the relationship between “weak” 
and “powerful” catalysts (possessing small or large values of rrii) leading to the GBP supercrit¬ 
icality, criticality or subcriticality. 

Clearly, C C {(mi,..., mAr) G : det{D — I) = 0} though the latter set also covers 
the case when A^ > 1, p{D) > 1 and D has an eigenvalue 1. Using the Laplace expansion of 
determinants we derive that for A^ > 1 one has det (D — J) = 0 iff 

N 

{airrii - 1) det{D - I)i^i = (1 - a*) ^ det {D - I)- - 

i=i 

where {D — I)^- is the {N — 1) x [N — 1) matrix that results from deleting the Uth row and 
the j-th column oi D — I, i,j = 1,... ,N. Since matrix D is irreducible by Lemma the 
maximal eigenvalue of any principal sub-matrix of D is strictly less than p{D) (see, e.g., [TU] . 
voL2, Ch.l3, Sec.3). Thus, if p{D) = 1 then det(Zi) — J)*,* ^ 0 for any i = 1,..., A^ and 

^ _ ^ ~ (^ ~ 

^ (1 - «,) (- 1 )-+^-V^F^„^^.(oo) det (D - 

aidet(D - I)i^i 

Obviously, the matrix (D(mi,... ,mjv) — does not depend on m* and rrij for any i, j = 
1,..., A^. Therefore, if (mi,..., rrin) G C and all rrij, for j = 1,..., N, j ^ i, are fixed then m* 
is uniquely determined by ([^. 

Let us hnd the bounds for possible values of m,, i = 1,..., N, such that (mi,..., m^) G C. 
At first recall that according to the Perron-Frobenius theorem, for irreducible matrices A and 
R, the Perron root p{A) does not exceed the Perron root p{B) whenever all the elements 
of matrix B — A are nonnegative. Furthermore, then A = B ii p{A) = p{B). Assume that 
m* = 0 for all i = 1,..., A^. In this case the maximal row sum of D is less than 1 and, 
consequently, p{D) < 1 in view of Theorem 1.5 and Corollary 1 in [26], Ch.l, Sec.l and 4. 
If now we let mi —)■ oo then the Perron root p{D{mi,0 ,... ,0)) strictly grows to infinity by 
virtue of the Perron-Frobenius theorem proof. Since p{D) is continuously dependent on any 
element of matrix D, there exists Mi > 0 such that p (R(Mi, 0,..., 0)) = 1. Moreover, if 
mi < Ml then p (Zl(mi, 0 ,..., 0 )) < 1 whereas for mi > Mi and any nonnegative m 2 ,..., mAr 
one can verify that p (Zl(mi, m 2 ,..., > 1. The exact value of Mi can be found by setting 

m 2 = m 3 = ... = mAr = 0in([^. In a similar way, we may indicate positive M 2 ,..., Mat 
such that (mi,..., rriN) G C with rrii = Mi and rrij = 0 for each i = 2,... ,N and j 7 ^ i. 
Thus, we have found the smallest parallelepiped [0, Mi] x ... x [0, Mat] containing C, i.e. 
C C [0,Mi] X ... X [0,MAr]. 
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Figure 1: Examples of the criticality set C for N = 2 and N = 3. 


Assume now that mi = m^ where G [0,Mi). Reasoning as above we conclude that 
there exists M 2 {m\) G (0, M2] such that the equality p M2(m'j’), 0,..., 0)) = 1 holds 

true. The value M 2 (m\) is given by formula Q when mi = m\ and m3 = ... = mjv = 0. 
Let m 2 = m^ with m^ G [0,M2(m5)). Next we can take M^{m\,m'^ G (0, M3] so that 
p (Zi)(m5, m^, M3(m3, m^), 0,..., 0))) = 1. Using identity Q and by way of the above rea¬ 
soning we choose the array m° G [0, Mj(m5,... ,m°_3)) and thus reveal a collection of numbers 
Mj(m5,... ^’mPi-i) e (0,Mj] such that p [D{m\,... ,m°_3, Mj(m‘j’,... ,m°_3),0 ,... ,0)) =1 for 
z = 1,..., — 1 (naturally, we set Mi(0) = Mi). Finally, to pick m^ for (m?,..., m^) G C we 

employ ([^ with mi = m\, ..., miv-i = ^%-i- Thereby we demonstrate how one can vary a col¬ 
lection (mi,... ,m^) providing the criticality of the corresponding CBP. Note that if for some 
step z = 1 ,..., — 1 we choose m° = Mi(m5 ,..., m°_i) then m° = 0 for all A; = z -|- 1 ,..., A^. 

Moreover, if for some step z = 1,..., A^ — 1 we take m^ > Mj(m5,... ,mi_i) then CBP is 
supercritical for any nonnegative values mj+i,..., m^r. At last, if m^ > 0 then the choice 
(mi,... ,m^_i,mAr) with mjq > m% or m^ < m% leads to supercritical or subcritical CBP, 


respectively. 

To illustrate this discussion we provide below Figure [T] depicting the set C for some particular 
cases. On the hrst picture there is a black curve representing the set C for a simple random 
walk on Z with two catalysts located at neighboring points. In this case, for instance, by 
Theorem 2 in P, we deduce that w = 1/2 for each i,j = 1,2. We also set ai = 0.3 

and a 2 = 0.8. The domain marked by light grey color under the curve is related to subcritical 
CBP whereas the domain marked by dark grey color over the curve corresponds to supercritical 
CBP. On the second picture there is a surface depicting the set C for a simple random walk 
on Z with three catalysts located at three subsequent points. With the help of Theorem 2 in 
[H] we derive that WjFwi,Wj{oo) equals either 1/2, given that \wi — Wj\ < 1, or 0, provided that 
\wi — Wj\ = 2, i,j = 1, 2,3, except for the case = 0- We also assume that ai = 0.3, 

Q!2 = 0.6 and = 0.8. The domain under the surface corresponds to subcritical CBP while 
that over the surface is related to supercritical CBP (the net on the surface displays the level 
lines). 
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4 Moment analysis of CBP 


It was Ch.J. Mode who showed (see [22]) that the mean particles numbers in a supercritical in¬ 
decomposable multi-type Bellman-Harris process grow exponentially (as time tends to inhnity) 
with a certain rate v called a Malthusian parameter. The parameter v is positive and deter¬ 
mined as the unique solution to equation p(H{u)) = 1. Lemma demonstrates that for our 
auxiliary L-dimensional supercritical Bellman-Harris process this procedure can be simplihed 
and the Malthusian parameter can be found as the unique solution to equation p{D{v)) = 1 
(observe that square matrix H{\) is of size L > N whereas square matrix D{X) is of size N). 
So, in the framework of CBP the Malthusian parameter u is dehned as the unique (positive) 
solution to equation pi^Diy)) = 1 whenever p{D) > 1. Given p{D) = p(Zi)(0)) > 1, such so¬ 
lution always exists, since according to the Perron-Frobenius theorem the function p{D{X)) in 
variable A is strictly decreasing and tends to 0 as A —)■ cxd. 

Recall that, for an n x n matrix R, the matrix Rij or {R)ij is the (n — 1) x (n — 1) matrix 
that results from deleting the Tth row and the j-th column of R, i,j = 1,, n. To formulate 
our main result we need some new notation. For supercritical CBP consider functions a„(a;, y), 
x,yES, ugN, which appear in description of the asymptotic behavior, as f —>■ oo, of the u-th 
factorial moments m„(f; x, y) := Ea;/x(f; y){p(t; y) — 1)... {pit] y) — nRl) of the local particles 
numbers in CBP (index x denotes the starting point of CBP). The function ai{x,y) has the 
most concise form when x,y G IF, namely. 


ai{wi,Wj) 




{v + i3j)^'{vy 


i,j = 1,...,N, 


( 8 ) 


where Ajj(A) := (—1)*+-^ det(/ — D{X))ij and A'(i/) := ^ det(/— iA(A)) If iV = 1 then 
All (A), A > 0, is supposed to be identically 1. For x E W and y E S \ W we set 


ai{wi,y) 


EjL,A,,(zz)(l-a,)q(zz)i^F;.,^(zz) 

{^ - Q{y,y))^'{A{^ - 


i = 1 ,..., A. 


(9) 


Here HFy y,{X) := {-q{x',x')){X-q{x',x')) VF*,y(A), A > 0, x',y' E S,H C S, i.e. //Fi,/y(t), 
f > 0, is the c.d.f. of a hitting time of state y' under taboo H given that the Markov chain r] 
starts at state x' (see, e.g., 0 )- Note that G* (A) = (3j{X+j3j) ^ for each A > 0 and j = 1 ,..., A. 
The values an{x,y) for n > 1, a; G IF and y E S are evaluated with the help of the following 
iterative scheme 


N 

k=l 


/3kAkAnA 

{nu + j3k)A{nv) 


hn,A^i{wk,y), ■ ■ . ,a„-i(wfc,2/)) 


( 10 ) 


where i 


1 ,..., A, A(A) := det(/ - D{X)), A > 0, 

FT) 




E 


n\ 


r=2 


r\ ^ 

2lH-...+2r=n 


( 11 ) 


and zi,..., Zn-i > 0, fc = 1,..., A. In this way the values a„(a;, ?/) for n > 1, x G S' \ IF and 
y E S are determined by formula 


aAx,y) = 2^.^^WiFy^Anu)an{wi,y). 


( 12 ) 
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Define functions An(x), x G S, n G N, arising in the asymptotic behavior of the n-th 
factorial moments M„(f; x) := Exfj,(t){fj,(t) — 1 )... — n + 1) of the total particles numbers 

in supercritical CBP. Namely, for n = 1 and x G W, put 


N 


Mwi) = 


z/(z/ + /3j)A'(z/) 


i = 1, 


.N. 


(13) 


When n > 1 and x & W the values An{x) are evaluated according to the following iterative 
scheme 


N 




An{w^) ^ 


^n,j^A\(Wj^^ . . . , (iCj)) 


(14) 


where i = 1,..., iV. Bearing on the latter equalities the values An{x) for n > 1 and x & S\W 
are obtained by way of 

^nix) = W,F*n,^{niy)AniWi). (15) 

For critical CBP introduce functions hn{x,y) and Bn{x), x,y E S, n E N, appearing in the 
asymptotic behavior of mn(t;x,y) and Mn(t;x), respectively. Namely, for i, j = 1,..., A and 
yES\W, set 


, / X Aj j(O) , , ^i,*(0)(l Oij)wFwj,y{oo) 

bi{wi,Wj) = bi{wi,y) = - - - - XX , XX 

/3jA'(0) (-g(l/,l/))A'(0)(l - wFy^yipo)) 


(16) 


where A'(0) := Adet(/ —D(A)) When n > 1, i = 1,..., N and y E S, the values 


bn{wi,y) are evaluated by the iterative scheme 


bn{wi,y) — ^ 2(n - l)A'(O) / ^ ]br{wk,y)bn-r{wk,y). 


(17) 


Thus the values bn{x, y) for n > 1, x E S \ W and y E S are determined by the formula 


E n 

. w,F^,wi{^)bn{wi,y). 

1=1 


(18) 


For f = 1,..., A put 


N 


Bi(wi) = 

1=1 


ai(/j(l) - 1 )A,,( 0 ) 


A'(0) 


When n > 1 the values Bn{wi) are computed iteratively, namely. 


(19) 


N 


a,/;(l)A,,< 0 ^^/n\ , jsr 

qni i r / I 1,..., A. 


Bniw.) Y 2(2„ - l)A'(O) ^ \r 


( 20 ) 


Bearing on the latter equalities the values Bn{x) for n > 1 and x E S \ W are given as follows 


E n 

WiFx Wi (^^).®n(^i) • 

1=1 


( 21 ) 
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At last, for subcritical CBP define functions Cn(x), x E S, n E N, arising in the description 
of the asymptotic behavior for Mn{t;x). Set 


N 


Ci{wi) = l + ^ -^ % = 


i=i 


A(0) 


.N 


.N 


Ci(a;) = l-y^ w,F^,wi{oo)+ y^ WiF^^^.{oo)Ci{wi), xeS\W, 
and, for n > 1, put 


N 


c„{w,) = Y, ■ ■ ■ ,Cn-l(n)), ! = 1,...,iV, 

7 = 1 ^ ^ 


( 22 ) 

(23) 

(24) 


E N 

W,F:^,Wi{oo)Cn{Wi), X E S \W. (25) 

1=1 

The following theorem is the main result of Section [^providing the moment analysis of the 
local and total particles numbers in CBP. This theorem generalizes Theorems 4.1, 4.2 in [3l] 
and Theorem 1 in [TT] as well as some statements of Theorem 2 in [2]. 


Theorem 1 For each n G N, given that < cx) for any i = 1,... ,N, functions mn{t] x, y) 
and Mn(t;x) are bounded on every finite interval in [0, oo) for fixed x,y E S. Moreover, 
whenever Eer < oo for any i = 1,... ,N, the asymptotic behavior of mn(t;x,y) and Mn(t;x) 
as t ^ oo depends essentially on the class of CBP and is as follows. 

1. If p{D) > 1 then z/ > 0 and 

mn{t-,x,y) = an(a:,|/)e”''* + o (e”"'*) , (26) 

Mnit;x) = Anix)e^‘'^ + o {e^^'^) , t ^ oo, (27) 

where functions an{x,y) and An{x), x,y E S, are strictly positive. 

2. If p{D) = 1 then 

mn{t-,x,y) = hn{x,y)C~^ + o{C~^) , (28) 

Mn{t;x) = Bn{x)f^"'~^ + o , t ^ oo, (29) 

where function bn{x, y) is strictly positive for any values x,y E S iff u dwjFwi,wj (u) < oo 
for each i,j = 1,...,N and, in addition, in case of recurrent p and n > 2 one has 
^^1 Oj |/j(s) — s| > 0 for some s E [0,1). Moreover, Bn{x) is strictly positive for 
any x E S iff the Markov chain p is transient and udw^Fy^.^^jiu) < oo for each 
i, j = 1,..., iV. Otherwise, respectively bnf, •) = 0 and Bnf) = 0. 

3. If p{D) < 1 then 


mn{t-,x,y) = o(l), (30) 

Mn{t;x) = C„(a;) + o(l), t-)■ oo, (31) 

where function Cn{x) is strictly positive for each x E S iff the Markov chain p is transient. 
Otherwise, C„(-) = 0. 
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Stress that we do not impose any restrictions on the character of the Markov chain rj such 
as symmetry and homogeneity of transition probabilities or hnite variance of jump sizes etc. 
Therefore, in Theorem we establish exact asymptotic relations in item 1 and only partially 
in items 2 and 3. It is not possible to give any useful and general information on items 2 and 3 
without further knowledge of the underlying motion. For instance, even for critical symmetric 
branching random walk on with a single catalyst (see 0 ) there are four different asymptotic 
formulae for m„(f;x,|/) and Mn(t;x) depending on dimension d = 1,2,3 or 4 and thereby 
on the decay rate of transition probabilities (note that for d > 5 Theorem generalizes the 
corresponding results in [2]). Moreover, for subcritical branching random walk on with a 
single catalyst (see [7]) the decay orders of mn{t]x,y) do not coincide for different d G N. 

Let us proceed to the proof of Theorem [1} 

Proof. We will mainly employ the auxiliary Bellman-Harris process constructed in Sec¬ 
tion Firstly we establish Theorem for x,y E W (Case 1) and afterwards extend these 
results to the general case x,y E S (Case 2 ). 

Case 1. So, assume that x,y E W. The following proof of Theorem[l]for Case 1 is divided 
into 9 steps. At step 1 we derive a system of renewal equations involving the factorial moments 
of the local particles numbers in CBP. Then at steps 2, 3 and 4 we employ this system to prove 
the corresponding assertions (26), (28) and (30) of Theorem related to moment analysis of 
the local particles numbers. Next at steps 5 and 6 we derive a system of renewal equations 
involving the factorial moments of the total particles number in CBP depending on whether 
the Markov chain rj is recurrent or transient. At last, at steps 7, 8 and 9 we exploit these 
systems to establish statements (27), (29) and (31) of Theorem 1, respectively, concerning the 
moment analysis of the total particles number. Note that at step 6 we introduce an auxiliary 
decomposable {L + l)-dimensional Bellman-Harris process with a final type of particles. 

Step 1. At this stage we derive a system of renewal equations involving the factorial mo¬ 
ments of local particles numbers in CBP. Denote by k, 1) := EikZi(t){Zi{t) — l)... {Zi(t) — 

n -|- 1), for f > 0, n G N and k,l = 1,..., L, the n-th factorial moment of the number of par¬ 
ticles Zi{t) of the l-th type at time t in the Bellman-Harris process (index k stands for the 
type of the parent particle). According to the construction of the auxiliary process one has 
rrinit] Wi, Wj) = i,j) for each n E N, t > 0, i,j = 1,..., N. Thus, given that < oo 

for some n G N and any k = 1,... ,N, the hniteness of functions mn{t] Wi, Wj) for each t > 0 
and i,j = follows from Theorem 1 in [27], Ch. 8 , Sec. 6 . Moreover, in view of this 

theorem functions k,l), t > 0, k,l = 1 ,..., L, satisfy the following system of renewal 

equations 


mf" (t; k, 1) = 6k,iil - Gfc(t)) + 


^ pt 

r=l Jo 


{t — u] r, 1) dGk{u) 


while for n > 1 this theorem and the generalized Faa di Bruno’s formula (see 

m^^{t — u; r, /) dGk{u) 


entail 






r.l •'0 


(32) 


( 33 ) 


n\ 






X 


1.0 

p\ 


ds 


J{L) 


s=(l,...,l) 



;( 1 ) 

■P 


rripit - U] L, 1) 
p\ 


tp 


dGkiu) 


12 













Here the symbol Y2 nieans the sum taken over all L-dimensional vectors ii = 


'1 ) 




irt-i = with nonnegative integer components satisfying the following equality 

Z]p=i P (^4^^ + ... + 4^^ j = n. For every r = 1,..., L, we also put J(r) := X]p=i 4^^- Substi¬ 
tuting the explicit formulae Q an d ([^ for the offspring generating functions gk{ ), k = 1,..., L, 
of the Bellman-Harris process in (|33| we come to more concise relation 


m. 


BH 


^ nt 

= ''S^mk^r / {t — u;r,l) dGk{u 

r=i ’ 


(34) 






p\ 


where !(•) is the indicator and the symbol 'Yh means the sum taken over all nonnegative integer 


H, • • in-i such that Y!l=iP‘^p — use the notation J := H + 


+ 'in-l- 


Step 2. Now we prove the part of assertion ( |26[ ) related to the case x,y & W. So, assume 
that p{D) > 1. Applying Theorem 2.1, item (iii), in [T2] to the system of renewal equations 


(32) we get 




(-1)^+1 det(/-iJ(z/)),-, 

(^ + /5,.)|.det(/-i7(A))L^ 


t —)■ oo. 


(35) 


\=u 


for each i, j = 1,..., A^. Here u > 0 and ^ det(/ — if(A))|^_^ > 0 (see, e.g., [22])- Comparing 
formulae (j^ and (26) with (35) we see that one needs to establish relations between determi¬ 
nants and algebraic adjuncts of matrices I — H{X) and I — D{X). The following lemma provides 
them. 

Lemma 2 For each A > 0 and i = 1,..., \Kj\, j, k = 1,..., N, one has 

det(J - H{X)) = det(J - D{X)), det(J - H{X))j^k = det(/ - T>(A))j-fc, 

(_l)HT+m det(/ - H{X))LU)+i,k = (-1 )^+Xto)+*G';(A) det(/ - D{X)),,k. 

Proof. We apply transformations to the columns of matrix I — H{X) which do not change 
its determinant. Namely, for each j, k = 1,..., N add the {L{j) -|- z)-th column multiplied by 
^L(j)+iW to the k-th column if k E Kj, i.e. k = k{i,j) for some i = 1,... ,\Kj\. After 
these transformations we get a block matrix consisting of four blocks. The left upper block 
of size At X At is just I — D{X) whereas the left lower (L — N) x N and the right lower 
(L — At) X (L — At) blocks are the zero and the identity matrices, respectively. Employing the 
formula for the determinant of a block matrix (see, e.g., [12], vol.l, Ch. 2 , Sec.5) we come to 
the hrst assertion of Lemma The second one is established in the same manner. To prove 
the third assertion we apply the above transformations to the matrix (/ — H(X))L{j)+i,k and 
then shift the column number L{j) + i — 1 to the column number N. The determinants of the 
matrix (/ — iL(A))LQ)+i,fc and the transformed one are equal up to the factor (—l)'^ 0 ')+*--'v-i_ 
Afterwards, we involve the formula for the determinant of a block matrix once again and use 
the Laplace expansion of determinants with respect to the W-th column. Lemma is proved 
completely. □ 

Thus, by virtue of (35) and Lemma we prove (26) when n = 1 and x,y E W. Before 
focusing on the case n> 1 let us note that in [22] Ch.J. Mode found the asymptotic behavior 
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of the mean particles numbers in supercritical multi-type Bellman-Harris process under rather 
restrictive conditions that each c.d.f. Gfc(-) has a square-integrable density. This is the reason 
why we bear on similar results by K. Crump established in [12] without any additional assump¬ 
tions and not on the mentioned results by Ch.J. Mode. To check relation (26) when n > 1 and 
X, 1 / G hh we employ the induction method in variable n. The case n = 1 is verihed above. Let 
formula (26) be valid for all the moment orders not exceeding n — 1. Then the second term at 
the right-hand side of (34) denoted by k, 1) has the following asymptotic behavior 


v;(«; ]V})e"-‘GMnO - 




n—1 


ll\ . 


n 

p=i 


p\ 


as f —)■ cx), for each k = 1,... ,L and I = 1,... ,N. Taking into account the latter relation and 
applying Theorem 2.1, item (iv), in [12] to the system of renewal equations (34) we obtain 




N 


-l)^+^/jfc <let{I-H{np))k,i 

^ (nz/ -f /dfc) det(/-iL(uz/)) 




as f —>■ cxD, for each i,j = 1,... ,N. Here we use the alternative representation (see, e.g., [T], 
Theorem 3.3) for function hn,k dehned in ( pT| . The latter asymptotic relation combined with 
Lemma 1^ and formula (10) leads to the desired statement in (26) when n > 1 and x,y E W. 
Note that functions an{wi,Wj), n > 1, i,j = 1,... ,N, are strictly positive according to, e.g., 
[H], vol. 2 , Ch.l3, Sec.3. 

Step 3. Next we prove assertion ( |2^ for x,y E W. Assume that p{D) = 1 or, equivalently, 
that the unique solution to p{D{v)) = 1 is z/ = 0. Applying Theorem 2 . 1 , item (iii), in [T2] to 
(32) we come to the following relation 

(l)*+Jdet(/-if( 0 )),, 




/3, |^det(/-Lr(A))|,^, 




.N. 


A=0-l- 


This formula along with (16) and Lemma entail the required statement of (28) when n = 1. 
For n > 1 we employ the induction method in variable n. Performing the induction step we get 


Vn{t;k,l) 


E iX 

i 


ockfk^iX)n\ fbp{wk,wi)y 


n— 1 


ip. . ..In-P 


n 


p=i 


pi 


{t - u)^-^dGk{u) 


n—1 


^n-2 ^2 

^ r=l 


n 


br{Wk,Wl)bn-r{Wk,Wl), t ^ OC, 


ioT k,l = 1,..., N, whereas Vnit; k,l) = 0, for k = N + 1,..., L, I = 1,..., N. Hence, applying 
Theorem 2 . 1 , item (v), in [T2] to (34) we see that 


i,j) A (- 1 ) ( 1 ) det (/ - if ( 0 )) ^ / rz 


kn—l 






0+ r=l 


y^i ]bJ^Wk,Wj)bn-JyWk,Wj) 


as f —?• oo, for each i,j = 1,...,N. The assertion (28) for n > 1 now follows from the 
latter relation, dehnition 0 and Lemma In view of [19], vol.2, Ch.l3, Sec.3, the algebraic 
adjunct Ajj(O) is strictly positive for any i,j = 1,..., N. As was established in [ 22 ] one has 
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A'(0) G (0, cxd]. Moreover, A'(0) 
Therefore, functions b„ 


< oo 


iff ud 




only if udwjF^.^^,^{u) < 


wi,wj{u) < oo for each z, j = 1 ,..., A. 
are strictly positive for any n G N, i,j = if and 

Wi,Wj\^j ^ oo and, in addition, in case of recurrent rj and n > 2 one has 
l/*('S) — s| > 0 for some s G [0,1). Note that the latter condition allows us to separate 
the case of ordinary motion of particles without breeding. Such CBP is critical for recurrent 
Markov chain rj and subcritical for transient rj. 

Step 4. Relation pOj ) for the case x,y E W is verihed by application of Theorem 2.2, item 
to (32) and (34). Here we essentially bear on the fact that M” converges elementwise 


(11 j, in _ _ 

to the zero matrix, as n —)■ oo, since in subcritical case p{M) < 1 on account of Lemma 

Step 5. Let us derive a system of renewal equations involving the factorial moments of 
the total particles number when x E W and rj is recurrent. Denote by Z{t) := J2j=i 
t > 0, the total particles number at time t in the auxiliary Bellman-Harris process. According 
to the construction of this process the distribution laws of p(t) and Z{t) are the same for 
each f > 0 iff the Markov chain rj is recurrent. Indeed, whenever the Markov chain rj is 
transient, with positive probability there are particles at time t in CBP having infinite life- 
lengths since they will never hit the set W after time t. These particles are not comprised 
by our L-dimensional Bellman-Harris process. So, we will now concentrate on the recurrent 
case. Then one has Mnit;Wi) = for each f > 0, n G N and i = 1,..., A where 

k) ;= EfcZ(f)(Z(f) — 1)... (Z(t) — n -(-1), f > 0, n G N and k = 1,... ,L. Thus, given 
that E^r < oo for some n G N and any Z = 1,..., A, Theorem 1 in |27], Ch. 8 , Sec. 6 , implies 
the finiteness of functions for each f > 0 and i = 1,..., A. Moreover, in view of 

this theorem functions fc), t > 0 , A; = 1 ,..., L, satisfy the following system of renewal 

equations 

^ M^^{t-u-,r)dGk{u). (36) 


Mf 


'(t; A;) = 1 - Gfc(A) + ^ mk,r f 

r=l Jo 


For n > 1 this theorem and the generalized Faa di Bruno’s formula (see PD!) entail equations 
in M^^{t-,k) the same as those obtained from (33) after replacing [t] k,l) by [t] k). 
Therefore, for k), t >0, k = 1,..., L, the system of renewal equations (34) holds true 

if function k, 1) is replaced by k). 

Step 6. Next consider CBP with the underlying motion of particles governed by transient 
Markov chain rj. To cover the transient case we construct a new auxiliary Bellman-Harris 
process with particles of L -|- 1 types such that the (L -|- l)-th type is final (see, e.g., |27], Ch.5, 
Sec.3, and ra)- For the new Bellman-Harris process, the c.d.f. of the life-length of a particle 
of j-th type is denoted by Gj{t), A > 0, whereas its offspring generating function is gj{s, sl+i), 
s G [0,1]^, sl +1 E [0,1], j = 1,..., L -|- 1. Let Gi{t) = Gi{t), t > 0, for z = 1,..., L, and 


gj{s,SL+i) = gjis), s G [0,1]-^, 
in ([^ and ([^. Then put 


sl +1 G [0,1], for j = a -f 1,..., L, where Gi and gj are defined 


N 


gji^y ^L+l) T (1 ^ ^ 


Wk ^ Wj ,Wk 


(O)sfe 


k=l 


\Kj\ f ^ \ 

-|-(1 — aj) 'y ^^ Ti^j{po)sL{j)+i -|- (1 — Qj) I 1 ~ y ' ^ , WkFwj,Wk{’^) j '5L+1) 

i=l \ k=l J 

for j = 1,...,A, s G [0,1]'^, sl +1 E [0,1], i.e. definitions of functions gj and 'gj differ by 
the last term only. Since the {L -|- l)-th type is final, that is every particle of the {L + l)-th 
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type has infinite life-length and does not produce any offsprings, for the sake of definiteness 
we may set G£,+i(f) = 0, f > 0, and ^l+i(s,sl+i) = s G [0,1]^, G [0,1]. The mean 
matrix M = process has entries fhk^i = rrik^i, for fc, / = 1,..., L, fhi^L+i = 

(1 - ai) (^1 - Y.k=i WkFwi,iOki^)j, i = ^, ■ ■ ■ ,N, TriL+i,L+i = 1 and fhk^i = 0 for the remaining 

pairs of k and 1. Let Zj{t),t > 0, j = 1,..., L -|- 1, be the number of particles of type j at time t 
in the new Bellman-Harris process. Evidently, the distributions of vectors {Zi{t),... ,Zi{t)) and 
{Zi{t),... ,ZL{t)) coincide for each f > 0 whenever the parent particles of both Bellman-Harris 
processes have the same type. Moreover, the new process takes into account even the particles 
from the {L + l)-th group of particles in CBP (see Section which are not comprised by the 
L-dimensional Bellman-Harris process. Thu^for each n G N, i = 1,..., N and t > 0, one 
has Mn{t;Wi) = where we write fc) := EfcZ(f)(Z(f) — 1 )... (Z(f) — n-|- 1 ), 

k = 1,..., L + 1, and set Z{t) := ^^^=1 Zj{t). Then applying Theorem 1 in [27], Ch. 8 , Sec. 6 , 
we ascertain that whenever < oo for some n G N and any / = 1,..., A^, functions M„(f; Wi) 
are finite for each f > 0 and z = 1,..., iV. This theorem also implies that functions k), 

t > 0, k = 1,..., L + 1, satisfy the following system of L -|- 1 renewal equations 

^ ^+1 pt _ ^ 

k) = 1- Gkjt) + y^^mk,r / - u]r)dGk{u). 

r=l 

The {L -|- l)-th equation is just M^^(t]L -|- 1) = 1, f > 0. Substituting that value 1 and 
expressions for rfik,r, r = l,...,L-(-l, into the first L equations we get a new system of L 
renewal equations 


_ ^ pt 

k) = rrik^r / M^^{t — u; r) dGk{u) 

r=l ’ >^0 


(37) 


N 


+ 1 — Gk{t) +l{k E {!,..., N}) (1 — offc) I 1 — WjFwk,wj{'^) ) Gk{t) 

i=i 

for f > 0 and fc = 1,... ,L. It is not difficult to verify that in view of Theorem 1 in 
Ch. 8 , Sec. 6 , the system of L -|- 1 equations resulting from (33) upon replacement of L by L -|- 1 
and of (t; k,l) by M^^{t]k) holds true for n > 1, f > 0 and k = 1,...,L -|- 1. The 
(L -I- l)-th equation is just L -|- 1) = 0. Substituting that in the first L equations and 

simplifying them similar to derivation of (34) we get the system of L renewal equations in 
functions k), for n > 1, f > 0 and k = 1,..., L, analogous to (34) but with /c, 1) 

replaced by k). 

Step 7. Now we prove assertion ( [^ for the case x,y E W. So, assume that p{D) > 1. At 
first consider the recurrent case. Applying Theorem 2 . 1 , item (iii), in [T2] to (36) we see that 

r \ 

(38) 


oVt 


N 


^z/A'(z/)' 
i=i ^ ^ 


N 


1 - ajG*{iy) - (1 - aj)G*{iy)J2w^F.., „.. [u 


Wj,Wk ' 


k=l 


AiJu) 


N N 

uA'(u) 

= 1 .7 = 1 ^ ^ 


Sj,k - ^j,fcai/)(l)G'*(z/) - (1 - <^j)G*{i')wkF^ 


k=l 3= 

"«,(/;(!) - i)G-(OA,,i(o A «iWm - 


E 

i=i 


uA'A) 


E 

j=i 


iy{iy + l3j)A'{iy) 
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as t —>■ oo and i = 1,..., Here we rely on Lemma and formulae 0-0 as well as on the 
Laplace expansion of determinants and obvious equality J2k=iWkFwj,Wki.'^) = 1 valid in the 
recurrent case. We also take into account that A(n) = 0 since 1 is an eigenvalue of matrix 
Relation (38) along with (13) lead to the desired statement (27) when n = 1 and the Markov 


chain rj is recurrent. The transient case is treated in the same way except for two differences. 
The first one consists in employment of equations system ( [^ instead of ( [^ . The second one 
is that for transient Markov chain rj the strict inequality X]a:=i Wk^wj,wk{^) < 1 holds true at 
least for some j G {1,... ,iV}. However, when deriving relation for similar to (38), 

the gap between the latter sum and 1 is compensated by means of the additional summand in 


(37). Thus, when n = 1 assertion (27) is established and the asymptotic behavior of Mi(f;tCj), 
i = 1,..., A^, in supercritical CBP does not depend on whether the Markov chain rj is recurrent 
or transient. Since for n > 1 the systems of equations in k) and k) respectively 


are of the same type as (34), the asymptotic behavior of Mn(t]Wi) in supercritical CBP is also 


independent of recurrence or transience of rj. Moreover, when n > 1 derivation of (27) for the 


total particles number almost literally repeats that of (30) for the local particles numbers and 
thus is omitted. Note that function Ai{wi) (and, consequently, An{wi) for n > 1) is strictly 
positive for each z = 1,..., A^ in view of positivity of the terms at the right-hand side of (38). 
Step 8 . Next proceed to the proof of statement (^^ for the case x,y E W. So, consider 


the case p{D) = 1. Applying Theorem 2.1, item (v), in [12] to (37) we get 


N 




i=i 


(1 aj)Aj^i(0) 

A'(0) 


N 




oo 


(39) 


k=l 


N N 


%7oT - (1 - (^j)wkFw,,wk{oo)) 

k=i j=i ^ ' 


N 


^ ’ ’ 
1=1 ^ ' 1=1 


A'(0) 


for z = 1,..., A^. Here we use Lemma as well as the Laplace expansion of determinants 


and equality A(0) = 0 valid by virtue of assumption p{D) = 1. Relation (39) combined with 


(19) implies the desired assertion in (29) when n = 1 and Markov chain p is transient. When 


Markov chain p is recurrent, the corresponding assertion follows from Theorem 2.1, item 


in na and observation that the non-integral term at the right-hand side in ( p36| tends to 0 , as 
t —)■ oo, whereas that in (37) converges to a positive limit at least for some k G {1,..., A}. 


For n > 1 derivation of (29) for the total particles number again reproduces that of (28) for 


the local particles numbers and so is omitted. Note that in view of (39) function Bi{wi) (and, 
consequently, Bn{wi) for rz > 1 ) is strictly positive for each z = 1 ,... , A iff A'(0) < oo and 
Yl,k=i Wk^wj,wk{^) < 1 at least for some j G {1,..., A}. 

Step 9. At last, let us establish ([3l| for the case x,y E W. Thus, assume that p{D) < 1. 
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Applying Theorem 2.2, item (ii), in [23] to (37), we see that 




A(0) 


j=l \ k=l 

^ ^ A-rot - 

A(0) ~ - (1 - 

k=l j=l ^ ^ 

^ a,(f'(l) - 1)A,,(0) ^ ^ «,(/'(!) - 1)A,,(0) ^ 

E -- -i+E --. <-00. 


(40) 


A(0) 


for each i = 1,..., N. Here we involve Lemma as well as the Laplace expansion of determi¬ 
nants and inequality A(0) > 0 valid on account of assumption p(D) < 1. Formulae ( 22 ) and 


(40) entail the desired statement (31) when n = 1 and Markov chain rj is transient. When 


Markov chain p is recurrent, the corresponding relation ensues from (36) and Theorem 2.2, 
item (ii), in [23|. When n > 1, this theorem also implies (31) due to equations system (34). 


Observe that in view of (40) function C'i(wj) (and, consequently, Cn{wi) for n > 1) is strictly 
positive for each i = 1,..., A iff J2k=i WuPwj,w,,{oo) < 1 for at least some j G {1,..., N}, i.e. 
Markov chain rj is transient. 

Thus, Theorem [ 1 ] is proved completely for the case x,y E W. 

Case 2. Now we assume that either x E S'\hF or y E S'\hF or x,y E S \ W. The 
main idea of the rest of the proof is as follows. If a: E S' \ hF and y E W we supplement 
the set of catalysts W with x. Vice versa, if x G IF and y E S \ IF let us enlarge the 
set of catalysts IF by y. If x G S' \ IF and y = x then we will add x to the set of cata¬ 
lysts IF. If both X and y are from S' \ IF and, moreover, x ^ y, we supplement the set 
of catalysts IF with both states x and y. Afterwards we may employ the already estab¬ 
lished results for CBP with A -|- 1 or A -|- 2 catalysts. So, set wat+i = x, IF(x) := IF U {x} 
and IFj(x) := IF(x) \ {wi}, z = 1,..., A -|- 1. Let D{x; A) = {dij{x; A))^^\ be matrix with 
di,j{x;\) := 5ij aif'{l)G*{X) + {1- ai)G*{\)w^{x)Fli,w,W, ^ > 0. Here a^+i = 0, /)^+i(l) = 0 
and G^+iit) := 1 — t > 0. Similarly put 10^+2 = 2 /, IF(x,z/) := IF(x) U {y} and 

IF(x, z/) := IF(x,z/) \ {wj}, z = l,...,A + 2. Consider matrix A (x, z/; A) = {dij{x,y] 
for A > 0, such that dij{x,y;X) := 6ij aifl{l)G*{X) + (1 - ai)G*{X)w^ix,y)Fl,^,wjW- Here 
«v +2 = 0, /^+ 2 ( 1 ) = 0 and GN+ 2 it) := 1 — t > 0. Denote by z/(x) and Z 2 (x,z/) the 

unique solutions to equations p (iA(x; z/(x))) = 1 and p (iA(x, z/; z/(x, z/))) = 1, respectively. The 
following Lemmaimplies that u = v{x) = z/(x,p). 

Lemma 3 For each A > 0, the Perron roots p{D{X)), p(iA(x;A)) and p{D{x,y; X)) are simul¬ 
taneously greater than 1 , equal to 1 or less than 1 . 


Proof. The arguments of Lemma resemble those of Lemma It is enough to establish 
the desired statement for matrices D{X) and iA(x; A) only since automatically the similar as¬ 
sertion is true for D{x] A) and D{x,y, A). According to the Perron-Frobenius theorem matrix 
D(x; A) has a strictly positive left eigenvector v(x) = (z;i(x),..., z;Ar+i(x))''' corresponding to 
the eigenvalue p(iA(x; A)), i.e. (iA(x; A)"'^ — p(iA(x; A))/) v(x) = 0. We will apply the equiva¬ 
lent transformations to the obtained system of equations, namely, for each z = 1,..., A add the 
(A -|- l)-th row multiplied by WiFwN+i,wiW the Fth row. Focusing on the hrst A equations 
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only, we deduce that 


(ui, ...,vn) {D{X) - p{D{x] A))/) (41) 

= uat+i (p(-D(a;; A)) - 1) • • •, • 

Here we use the following identity 

Wi{x)P wj (•^) 

where i = 1,..., iV and j = 1,..., + 1. This equality is true, e.g., in view of Theorem 8 in 

|llj . Ch.2, Sec.11. It follows from (41) that 


(ui,... ,nAr) T)(A) > p(T)(a;; A)) (ui,... ,nAr) for p{D{x;X))>l, (43) 

(ni,...,nAr)T>(A) <p(T>(a;;A))(ni,...,n 7 v) for p(T)(a;; A)) < 1. (44) 


Examining the proof of the Perron-Frobenius theorem we conclude that due to (43) one has 


p{D{X)) > p{D{x]X)) for p{D{x-,X)) > 1. On account of Theorem 1.6 in [26], Ch.l, Sec.4, 
relation dil] ) yields that p{D(X)) < p{D{x;X)) for p(D{x;X)) < 1. These estimates entail the 
desired assertion of Lemma 1 □ 

After application of the established part of Theorem [T] to the CBP with the catalysts set 
W{x) or W{x^y) we observe that to complete the proof of Theorem one could bring the 
expressions like a„(r(;Ar+i, WAf+ 2 ) into the form (12). This can be easily realized by using the 
following result. 


Lemma 4 For x,y E S \ W, x y, i,j = 1,..., N and A > 0, one has 

det (/ - D(x-, A)) = A(A) (1 - wFlJX )), 
det (/ - D(x,A)) = A(A) (1 - „.F,y(A)) (1 - w(.)F;J\)) , 


— det(/-D(x;A)) 


— (l — 


x^x \ 


))> 


X=u 


J^det {I - D{x,y;X)) 


X=u 


(-1)‘+^ det (/ - D(x; A))^ . = A,,(A) (1 - ,„F;,,(A)) , 
(_l)j+-v+idet(/-£>(x;A))^„^, 




N 


1 - ivF;^^(A) 

(-l)‘+'*+idet(/-F>(x;A)) 




i=l 


N 


AT+I, —^ Aj'i(A)(l-aj)G*(A)wE (A), 

i=i 


det(J —iA(a:; A)) 


w+i.w+i 


N 


A(A) +5^iv.F;„(A)A,,,(A)(l-aj)G-(A)„,F;,JA), 

ti=i 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 
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(53) 


det(/ — D{x, y] X))n+ 2 ,n+i (l — wFy y{X)) 
(l - wF^^iX)) (l - w{x)Fly{X)) 


N 


= A{X)wF:^y{X) 


w,F* 


JA)A,,(A)(l-a,)G;(A)H^F;.,,(A) 


*j=i 


where v is the unique solution to pi^Diyo)) = 1 existing whenever p{D) >1. //z/ = 0 then 
^ det (/ — D{x; A)) means the right derivative at 0. 

Proof. Apply some transformations to the columns of matrix I — D{x] A) which do not 
change its determinant. Namely, for each i = 1,..., N add the {N + l)-th column multiplied 


by WiF^^^^ .uj^{X) to the Ath column. Afterwards using (|4^ and the Laplace expansion of 


the determinants with respect to the {N + l)-th row we come to (45). Relation (46) can be 


ioT X > u and A(z/) = 0. Equality (49) is verihed in the same way as (45). 


considered as ensuing from (45). Formulae (45) and (46) imply (47) and (48) since A(A) > 0 


Let us check the validity of (50). Obviously, the following identity is true 


Wkix)Fl^^{X) 


w 




- A > 0, 


(54) 


Apply some transformations to the rows of matrix {I — D{x] X))^ which do not change 
its determinant. Namely, for each i = 1,..., N — 1, add the A-th row multiplied by (1 — 
ai)G*{X)wF*^.^^{X)/ (l — wF*^^{X)) to the f-th row. Afterwards employing (42) , ([^ and the 
Laplace expansion of the determinants with respect to the A-th row we derive (|50p. 


Implementing the same transformations as for proving (45) and using the Laplace expansion 


of the determinants with respect to the A-th column lead to identity (51). Formula (52) is 


verihed by involving relation (42), the linearity property of a determinant with respect to rows 


and the Laplace expansion of the determinants. 


Relation (53) is established by the combination of the previous transformations which 


do not change the determinant of (/ — D{x,y, X))n^ 2 ,n+i- Firstly, similar to proving (45), 


for each i = 1,...,A, add the (A + l)-th column multiplied by Wi(x)F*j^_^^^^.{X) to the i- 
th column and take into account a counterpart of (42). Secondly, by analogy with veri¬ 


fying (50), for every i = 1,...,A, add to the z-th row the (A -|- l)-th row multiplied by 
(1 — ai)G*{X)wF^.^^{X)/ (l — ivA*a,(A)). Thirdly, apply the Laplace expansion of the determi¬ 
nants with respect to the (A -|- l)-th column and the (A -|- l)-th row. At last, employ formulae 


(42) and (54) as well as the following equalities 


w{x)F ^,,y(A) -h 


Ww(x)F* y{X) 


1 - wF^^iX) 

W(x)Fly{X) (l - wFy^yjX)) 

(l - wF*^^{X)) (l - wix)F*^y{X)) 


= wF 




(A)- 


1 - 


W(x)-‘ 




1 - WF*^y{X) 

= wF*y{X), A > 0, z = l,...,A, 


valid according to the proof of Theorem 2 in |9]. In result we come to (53). This completes the 
proof of Lemma □ 

Thus, using Lem ma we derive relations ( [^, (12), (15), ([Ts] ), (21), (23), (25) from 
(10), (13), (14), ( [l6| ), ([T^, (19), (20), (22) and^4p. Theorem [ij is proved completely. □ 
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5 Applications 

In section by virtue of Theorem [T] we established that the Malthusian parameter v plays a 
crucial role in the asymptotic behavior of both total and local particles numbers in supercritical 
CBP. Before considering some particular examples in the present section let us discuss the 

alternative way of evaluating v. For this purpose introduce matrices Z)(A) = (di 


N 


and 


N 




D{X) = idij{X) ] with the corresponding entries 

/ i,j=l 


d,,(A) := (aJ'{l)G*{X)-l + {l-ai)G:{X) 


dij(A) . 


A - q{wi,Wi) 

-q{wi,Wi) 

6,, K/'(1)G:(A) - 1) (1 - , A > 0, 

A - q{wi,Wi) 


«,/'(i)g:(A) 
1 


1 + (1 — aj)G*(A)- 


-q{wi,Wi 


.(A) 


Gx{wi,Wj) 


ai 




’^q{wi,Wi 


G*{X), A>0. 


Here G\{x,y) := e~^^p{t;x,y) dt, A > 0, is the Laplace transform of the transition prob¬ 

ability p{t;x,y), x,y E S, t > 0, of the Markov chain p. Note that function Go{x,y) = 
limA-i-o-i- Gx{x, y) is called Green’s function and is finite iff the Markov chain p is transient (see, 
e.g., Theorem 4 and Corollary 2 in m, Ch.2, Sec.10). Recall also that according to Theorem 3 
and relation (4) in [TT], Ch. 2 , Sec. 12 , one has 


Wi,Wj 


(A) = 


Gx{wi,Wj) 

Gxiwj^WjY 


W^,Wi 


.:(A) = 1- 


(A - q{wi,Wi))Gx{wi,Wiy 


(55) 


for any z, j = 1 ,..., A^, i 7 ^ j and A > 0. 

It follows from Lemmaand dehnition of the Malthusian parameter in section]^ that v can 
be found as the maximal A satisfying the relation det(Zl(A) — I) = 0. The lemma below permits 
us to hnd V as the maximal A being the solution to equation detZl(A) = 0 or the equivalent 
equation det D{X) = 0. Observe that the entries of matrix D{X) are expressed via the Laplace 
transforms of hitting times under taboo sets Wi, i = 1,... ,N, whereas the elements of D{X) 
are represented in terms of the Laplace transforms of hitting times without taboo. At last, the 
entries of D[X) involve the Laplace transforms of transition probabilities only. 

Lemma 5 For any A > 0 one has det{D{X)—I) = 0 if and only if det D{X) = 0 or, equivalently, 
when detiA(A) = 0. Moreover, these relations are true even for A = 0 whenever p is transient. 

Proof. Introduce matrices R{X) = {Ti,j{X))f.^.^ and R{X) = {Y,j{X))^-^.^, A > 0, with the 
corresponding entries riq{X) := 5ij + (1 — 5ij) F^^wri,j{X) := Gx.{wi,Wj). By formula 


(55), the j-th column, j = 1,..., A^, of matrix R{X) is obtained from the j-th column of R{X) 
by multiplying it with Gx{wj,Wj) where Gx{wi,Wj) > 0 in view of irreducibility of p. Hence, 
det i?(A) = 0 iff detR(A) = 0, for any A > 0. Let us show that deti?(A) Y 0 whenever A > 0. 
Similarly to the arguments of Theorem 8 in HU, Ch. 2 , Sec. 11, one can derive that 


Wi ,Wj 


(t) — WjRwi,wAt) + 


N 

E 

k=l,k^j ' 


WkRwi,Wkit u) dF^^,^^i^Yu), 


(56) 
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for each = Note that the latter relation has a natural interpretation. Namely, the 

path of the Markov chain from Wi to Wj can either avoid the set Wj or can hit some site from 
Wj and afterwards reach wj. Applying the Laplace-Stieltjes transform to (56) we get 


N 


k=l, k^j 


(57) 


Multiplying each side of equality (57) by Gx{wj,Wj) and taking into account (55) we see that 

S- ^ 

Gx{wi,Wj) - _ =^WkFl^,^^{\)Gx{wk,Wj), i,j = 

A qyWi^ Wi) 

Rewrite these identities in the following matrix form 


R{X) - T(A) = U{X)R{X) or, equivalently, R{X) (/ - U{X)) = T(A) 


where 


r(A) := (5y (A - q(wu . U(\) := (* 


j Wi,Wj 


(A) 


N 


*J=1 


A > 0. 


Whence we deduce that deti?(A) ^ 0 for each A > 0, since detT(A) is strictly positive. 

It follows from (57) that {D{X) — I) R{X) = D{X). Considering the determinants of the 
matrices at the left- 
assertion of Lemma 


land and the right-hand sides of the latter equality we come to the hrst 
Its second assertion is implied by representation dij{X)Gx{wj,Wj) = 


djj(A) valid for each A > 0 and i, j = 1,..., 77 by virtue of (55). 


All the above reasoning holds true even for A = 0 whenever r] is transient. Thus, Lemma 
is proved completely. □ 

Let us consider some applications of our results to the models studied earlier by different 
researchers. 

Example 1. Focus on a catalytic branching random walk on Z*^, d G N, proposed in [31]. 
This model is a particular case of CBP if we set r; to be a symmetric and space-homogeneous 
random walk on with a finite variance of jump sizes as well as put A^ = 1, tci = 0 G Z'^, 
cxi = a, /3i = 1 and /i(-) = /(■). Then we deduce the same criticality condition as used in 
[28] . i.e. H = 1 or, equivalently, af{l) -|- (1 — a) (l -|- g(0, ^(0, 0)) = 1. For recurrent 

Markov chain p, Gq^( 0,0) is assumed to be 0. Applying Theorem to catalytic branching 
random walk on Z*^ we come to Theorem 4.1 and Theorem 4.2 in [3l| as well as some statements 
from Theorem 5 in [2H] and Theorem 1 in [H|. 

Example 2. Consider catalytic branching process with a single catalyst (located, say, at 
some site w & S) studied in [Tl]. Here the underlying motion of particles is governed by an 
irreducible Markov chain p. Thus, such setting is less restrictive than that in Example 1. As 
shown in [T3|, the asymptotic behavior of total and local particles numbers is determined by the 
mean offspring number, produced by a particle at the presence of the catalyst, being less than, 
equal to or greater than Gq^{w,w) -|- 1, and transience/recurrence of p. We come to the same 
classification letting N = 1, wi = w, ai = {1 — q{w,w))~^, /3i = 1 — q{w,w) and /((I) = m 


in our CBP. Then, in view of (55) and evident formula Fiij^yj{oo) = limA-i.o+^(A), the value 
sign{D — 1) coincides with sign{m — Gq^{w,w) — 1). Stress that in contrast to [H] we do not 
assume the existence of all moments of the offspring number, i.e. finiteness of /{^^(l) for any 
A; G N. Some assertions of Theorem [T] in the case of CBP with a single catalyst are stronger 
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than the corresponding statements of Theorem 1 in |T3] (for instance, cf. relations (30) and 
(31) in onr Theorem and point iii)a) in Theorem 1 in [Hj). However, some statements of 
Theorem 1 in [H] are not covered by onr Theorem [^becanse they involve asymptotic estimates 
for the moments of particles nnmbers in terms of local times at level w of the Markov chain rj. 
The anthors of na do not discnss the asymptotic behavior of those local times. 

Example 3. Concentrate on the branching random walk on d G N, with several sonrces 
investigated in [22] • This model is a particular case of CBP such that is a symmetric and 
space-homogeneous random walk on with a hnite variance of jump sizes and the symmetry 
of the random walk fails only at a hnite set of points of Z'^. As established in [23], the rate 
of exponential growth of the particles numbers is the maximal positive A satisfying equation 
det D{X) = 0 (this agrees with our Lemma [^. However, the necessary and sufficient conditions 
of existence of this positive solution A were not found as noted in Conclusion in [22] • It is 
worthwhile to remark that such necessary and sufficient condition is provided in our present 
paper and by virtue of Lemma and Theorem it is just p{D) > 1, i.e. in our terms amounts 
to a supercritical regime of CBP. 

Concluding the paper, we would like to observe that our approach of combination of hitting 
times under taboo and auxiliary multi-type Bellman-Harris processes permits us to obtain 
and justify the effective classihcation of catalytic branching processes with multiple catalysts. 
The results of this paper are valid under minimal restrictions on the character of motion and 
breeding of particles. Thus they generalize the previous works on closely related subject. The 
developed approach can be employed in the subsequent study of other characteristics of catalytic 
branching process. 
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